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Abstract 

Dijet differential cross sections for the reaction e'^p e~^ + jet + jet + X in 
the photoprodnction regime have been measured with the ZEUS detector at HERA 
using an integrated luminosity of 42.7 pb“^. The cross sections are given for photon- 
proton centre-of-mass energies in the range 134 < W < 277 GeV. The differential 
cross sections as a function of the dijet mass, and of the dijet angular variables 
have been measured for 47 < < 160 GeV and compared to next-to-leading-order 

QGD calculations. The dijet events in the region 75 < < 100 GeV have been 

used to derive a 95% G.L. upper limit on the cross section for photoprodnction of 
(7e+p^e+zox < 5-9 pb. Upper limits on the photoprodnction of new heavy resonances 
decaying into two jets are also presented for masses in the range between 60 GeV 
and 155 GeV. 
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1 Introduction 


In hadronic interactions, the distributions of the dijet mass, and of the angle, 
between the jet-jet axis and the beam direction in the dijet centre-of-mass system, provide 
a test of QCD. In addition, they are sensitive to the presence of new particles that decay 
into two jets. At HERA, these tests of QCD and searches for new particles can be made in 
photoproduction via ep scattering at ~ 0, where is the virtuality of the exchanged 
photon. Two types of QCD processes contribute to jet production at leading order (LO) 
in photoproduction |l|,||: either the photon interacts directly with a parton in the proton 
(the direct process) or the photon acts as a source of partons, one of which interacts with 
a parton in the proton (the resolved process). Jet photoproduction thereby provides a 
means to study the dynamics of both photon-parton and parton-parton interactions. 

At high Mjj values, the theoretical uncertainties due to hadronisation, multipartonic 
interactions and the limited knowledge of the photon parton densities are reduced. The 
dynamics of dijet production for hadron-induced processes has been investigated in detail 
in pp collisions at centre-of-mass energies of ^/s = 630 GeV and y/s = 1800 GeV 0. 
Next-to-leading-order (NLO) QGD predictions are in good agreement with the measured 
dijet mass and angular distributions up to 1000 GeV. Previous measurements 

of the dijet mass and angular distributions made use of cone algorithms for the 

identihcation of the jets. The kx cluster algorithm [|^|[ is used here, which allows a 
direct application of the theoretical jet algorithm to the data. The selection of high 
values permits a precise test of the description of the dynamics of dijet photoproduction 
to smaller distances than previously studied in photoproduction at HERA 0. 

Heavy particles that decay into two jets would show up as an enhancement in the dijet 
mass distribution. In particular, this mass distribution is sensitive to the production 
of the electroweak gauge bosons lE^ and Z^. New heavy particles decaying into two 
jets may also be identihed by deviations from the predictions of QGD in the Icos^*"^! 
distribution; in QGD, the | cos^*"^! distribution is peaked at unity, whereas many sources 
of new physics produce more isotropic angular distributions. 

In this paper, measurements of the dijet differential cross-sections da /and l/a da/dx, 
where y = (1 -|- | cos6'*|)/(l — | cos0*|), are presented as a function of M^K The variable 
9* is dehned by using the pseudorapiditie^ of the two jets with highest transverse energy 
in the event, cos^* = tanh|(77p*^ — and coincides with 9^^ for the case of 2 —> 2 

massless parton scattering. To increase the sensitivity to direct processes, measurements 
have also been made for dijet events in which the fraction of the photon’s momentum 
participating in the production of the dijet system is greater than 0.75. The NLO QGD 

^ The pseudorapidity is defined as ry = — In (tan |), where the polar angle, 6 , is measured with respect 
to the proton beam direction. 
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calculations ng are compared to the present measurements and upper limits are set 
on the inclusive cross sections for the photoproduction of and bosons. A search 
has also been carried out for the photoproduction of new heavy resonances decaying into 
two jets. 


2 Perturbative QCD calculations 


Perturbative QCD (pQCD) calculations of dijet cross sections in photoproduction can 
be written as a convolution of the subprocess cross section with the parton distribution 
functions (PDFs) of the photon and proton: 

dUep—>e jet jetX E I dXjf I dXpfh/pi^Xp, fJ^ppj d&ab—fjet jet(h'ij)) 

a,b do Jo Jo 


where y, x^ and Xp are the longitudinal momentum fractions of the quasi-real photon 
emitted by the positron, the parton a in the photon and the parton b in the proton, 
respectively; f^/e is the flux of photons in the positron and is usually estimated with the 
Weizsacker-Williams approximation |^; fa/^ ifb/p) represents the PDF of parton a {b) 
in the photon (proton) — in the case of direct processes, a is a 7 and fa/'y{x^, jjJp^) is 
given by 5(1 —x.^); the factorisation scale in the photon (proton) is denoted by yp^ i^Fp)] 
fiR represents the renormalisation scale; and the subprocess cross section, ddab^iet jet, 
describes the short-distance structure of the interaction. 


A wealth of data from hxed-target and collider experiments has made possible an accurate 
determination of the proton PDFs. In the case of the photon, experimental information 
on the quark densities is available from measurements of F 2 in e+e" collisions, while the 
gluon density remains poorly constrained. At high values of the dijet mass, the direct 
processes dominate and thus the sensitivity to the photon PDFs is reduced. 


The subprocess cross section is calculable in pQCD at each order. Recently, NLO QCD 
calculations in photoproduction ||9Hg have become available. In NLO QCD, the de¬ 
pendence of the calculations on the renormalisation and factorisation scales is reduced 
compared to LO. The results from the different NLO calculations have been compared 
and found to be in agreement within 5% in most of the phase-space region studied M . 


3 Experimental conditions 


During 1995-1997, HERA operated with protons of energy Ep = 820 GeV and positrons of 
energy = 27.5 GeV. The data sample used in this analysis corresponds to an integrated 
luminosity of 42.7±0.7 pb“^. The ZEUS detector is described in detail elsewhere [pA|p± 
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The main subdetectors used in the present analysis are the central tracking detector 
(CTD) positioned in a 1.43 T solenoidal magnetic held and the uranium-scintillator sam¬ 
pling calorimeter (CAL). 


Tracking information is provided by the CTD |]^ 


which is used to reconstruct the mo- 
< 9 < 164°. The relative transverse 


menta of tracks in the polar-angle^ region 15 
momentum, pt, resolution for full-length tracks can be parameterised as u{pt)/pt = 
0.0058 Pt © 0.0065 © O.OOld/pT^, with pt in GeV. The tracking system was used to es¬ 
tablish an interaction vertex with a typical resolution along (transverse to) the beam 
direction of 0.4 (0.1) cm and to cross-check the energy scale of the CAL. 

The CAL [|^ covers 99.7% of the total solid angle. The smallest subdivision of the CAL 
is called a cell. Energy deposits in the CAL were used in the jet Ending and to measure jet 
energies. Under test-beam conditions, the CAL single-particle relative energy resolutions 
were 0.18/(GeV) for electrons and 0.35 /E (GeV) for hadrons. Jet energies were 
corrected (see Section |j) for the energy lost in inactive material in front of the GAL. 

The luminosity was measured using the Bethe-Heitler reaction e~^p —*• e'^'jp 


19 . The 


resulting small-angle energetic photons were measured by the luminosity monitor, a lead- 
scintillator calorimeter placed in the HERA tunnel at Z = —107 m. 


4 Data selection and jet search 


A three-level trigger was used to select events online [1^,0. Events from collisions 


between quasi-real photons and protons were selected offline using criteria similar to 
those described in an earlier ZEUS publication [^. The main steps are briehy dis¬ 
cussed here. After requiring a reconstructed event vertex consistent with the nominal 
interaction position and cuts based on the tracking information, the contamination from 
beam-gas interactions, cosmic-ray showers and beam-halo muons was negligible. Gharged 
current deep inelastic scattering (DIS) events were rejected by requiring the total miss¬ 
ing transverse momentum, ^t, to be small compared to the total transverse energy, 

^t/ a/ Etf^ < 2 a/G eV. Neutral current DIS events with an identihed scattered-positron 
candidate [^| in the GAL were removed from the sample using the method described in an 
earlier publication The remaining background from neutral current DIS events was 
estimated by Monte Garlo techniques to be below 0.5%. The selected sample consisted 
of events from e'^p interactions with < 1 GeV^ and a median ^ 10“^ GeV^. The 


events were restricted to yp centre-of-mass energies in the range 134 < W < 277 GeV, as 
described in Section |^. 


^ The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the 
proton beam direction, referred to as the “forward direction”, and the X axis pointing left towards 
the centre of HERA. The coordinate origin is at the nominal interaction point. 
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The krp cluster algorithm was used in the longitudinally invariant inclusive mode [^] to 
reconstruct jets in the hadronic final state from the energy deposits in the CAL cells. The 
jet search was performed in the pseudorapidity (r/j-azimuth {(p) plane of the laboratory 
frame. The jet variables were defined according to the Snowmass convention |^. The 
jets reconstructed from the CAL cell energies are called calorimetric jets and the variables 
associated with them are denoted by V^caC There were 64708 events selected 

with at least two jets satisfying > 10 GeV and —1 < < 2.5. 


5 Monte Carlo simulation 


The Monte Carlo (MC) programs PYTHIA 5.7 [^] and HERWIG 5.9 were used to 
generate photoproduction events for resolved and direct processes. Events were generated 
using GRV-HO for the photon PDFs and MRS A for the proton PDFs. To 

study the dependence of the acceptance corrections on the choice of photon PDFs, the 
LAGl parameterisations were used. In both generators, the partonic processes are 
simulated using LO matrix elements, with the inclusion of initial- and final-state parton 
showers. Fragmentation into hadrons is performed using the LUND |M| string model as 
implemented in JETSET in the case of PYTHIA, and the cluster model in the case 
of HERWIG. For the measurements presented in this paper, the events generated using 
the PYTHIA and HERWIG programs were used for calculating jet-energy corrections and 
correcting for detector and acceptance effects. The corrections provided by PYTHIA were 
used as default values and those given by HERWIG were used to estimate the systematic 
uncertainties coming from the treatment of the parton shower and hadronisation. 


Photoproduction of the electroweak gauge bosons lU^ and was simulated with PYTHIA 
using the lowest-order processes qq and qq' —> lU^. Events were generated, includ¬ 

ing initial- and final-state parton showers, using GRV-HO for the photon and MRSA for 
the proton PDFs. 


To model the photoproduction of narrow heavy resonances decaying to dijets, MG events 
were simulated with PYTHIA using the LO process qq Z' and the same photon and 
proton PDFs as above. Samples of events were generated for each value of the mass from 
60 GeV up to 155 GeV in 5 GeV intervals. The vector and axial couplings to quarks 
and leptons were set equal to those of the Standard Model resulting in a width that 
increases linearly with mass. The cut | cos0*| < 0.6 was used to minimise the sensitivity 
of the acceptance corrections to the a priori unknown decay angular distribution of a new 
narrow state. 


All generated events were passed through the ZEUS detector and trigger simulation pro¬ 
grams based on GEANT 3.13 0. They were reconstructed and analysed by the same 
program chain as the data. The jet search was performed using the energy measured in 
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the CAL cells in the same way as for the data. The same jet algorithm was also applied 
to the hnal-state particles; the jets found in this way are referred to as hadronic jets. 


6 Jet energy correction and selection 

6.1 Jet hadron-level corrections 

The comparison of the reconstructed jet variables for the hadronic and the calorimetric 
jets in simulated events showed that no correction was needed for and 93 ^®* ~ 

and (pj®* ~ ‘^cai)- However, the transverse energy of the calorimetric jet underestimated 
that of the corresponding hadronic jet by an average of ~ 15%, with an r.m.s. of ~ 10%. 
The transverse energy corrections to calorimetric jets, averaged over the azimuthal angle, 
were hrst determined using the MC events. These corrections were constructed as factors, 
c(4T] which, when multiplied by provide the corrected transverse energies 

of the jets, 

6.2 Jet energy-scale corrections 

Further corrections to the jet transverse energy were developed to account for differences 
in the jet energy scale between data and simulations. This procedure relies on good 
understanding of the performance of the CTD track reconstruction in the selected region 
(see below). 

The response of the CAL to jets was investigated by using the following procedure pT| , ^ . 
In the central region, l^^j®*! < 1, the multiplicity distribution and the pt spectrum of 
charged particles associated with the calorimetric jets were compared for data and MC 
samples using the reconstructed tracks. The tracks were required to be in the ranges 
l^trackj ^ ^ > 300 MeV, where is the transverse momentum of the 

track with respect to the beam axis and is the track pseudorapidity. Tracks were 
associated with a calorimetric jet when the extrapolated track trajectory reached the CAL 
within a cone of one unit radius in the p-ip plane concentric with the calorimetric-jet axis. 
PYTHIA gives a good description of all the measured distributions. In this r/j®* region, the 
momenta of the tracks in the calorimetric jet were used to determine the total transverse 
energy carried by the charged particles, A'tL ark s) assuming zero mass for all tracks. Then, 
the ratio rtracks = -^TtracksZ-^Tcai Calculated and the distributions of this ratio for the 
dijet sample in data and MC generated events were compared. The mean value of the 
distribution in rtracks was determined as a function of pj®* for data (<'rtracks>data) and MC 
events (<'rtracks>Mc)- Differences between data and MC simulation of less than 1% were 
observed from the examination of the quantity (< rtracks > data / <rtracks>Mc) ~ 1- The 
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transverse energies of the calorimetric jets in the data were then modihed as a function 
of to correct for these differences. 

In the forward region, 1 < < 2.5, the energy scale of the jets was studied using 

the transverse-energy imbalance in dijet events with one jet in the central region and 
the other in the forward region. The distributions of the ratio raijet = 
jet)/i?;^* 3^1 (central jet) in data and the MC sample were compared. Differences between 
data and MC simulation of less than 2 % were observed from the examination of the quan¬ 
tity (<rdijet>data / <'rdijet>Mc) “ 1- The trausverse energies of the forward calorimetric 
jets in the data were then modihed as a function of 77 ^®* to correct for these differences. 
The widths of the distributions for rtracks and rayet are reasonably well described by the 
PYTHIA MC simulation, giving conhdence that the resolution in the energy of the jets 
is correctly described. 


The accuracy of the jet energy corrections was investigated by using data and simula¬ 
tions of neutral current DIS events at large where the transverse energy of the jet 
or jets is balanced by that of the scattered positron. The uncertainty in the absolute 
CAL energy calibration for the scattered positrons was 1% for positron energies above 
15 GeV [^. Jets were reconstructed in the laboratory frame using the algorithm de¬ 
scribed in Section The jet energy corrections described above were applied to both 
data and MC events. Then, the ratio r^/^et = where is the positron trans¬ 

verse energy, was formed and the distributions of this ratio for the inclusive jet sample 
in DIS data and MC simulation were compared. The mean value of the distribution of 
’"e/jet was determined as a function of r^j®* in the range —1 < r^j®* < 2.5 and of in the 
range 14 < E^^^ < 90 GeV for data (<re/jet>data) and MC events (<Te/jet>Mc)- Inspection 
of the quantity |(<re/jet>data /<'re/jet>Mc) — 1 | showed that the differences between data 
and MC simulation were smaller than 1%. This variation was therefore included as a 
systematic uncertainty in the present analysis. 


After these further corrections to the jet transverse energy, events with at least two jets 
satisfying E^,^^ > 14 GeV and —1 < r^j®* < 2.5 were used to measure the dijet cross sections 
presented in Section 


7 Reconstruction of jet and kinematic variables 

The invariant mass of the two jets with highest E^,^^ in the event was reconstructed using 
the formula 

Mjj = y^2E^®*^E^®*^[cOsh - 77^*2) _ pQg 

Only the absolute value of cos 6* can be determined because the originating parton cannot 
be identihed. For > 47 GeV and | cos0*1 < 0.8, the average relative resolutions in 
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Mjj and x were 8% and 7%, respectively; the average resolution in | cos 6*^ 


was 0 . 02 . 


The Jacquet-Blondel method |^, applied to photoproduction events was used to 
estimate W from the energies measured in the CAL cells, hhcai- Due to the energy lost in 
the inactive material in front of the CAL and to particles lost in the rear beampipe, WAai 
underestimates hh by ~ 10%, with an r.m.s. of ~ 5% This effect was corrected for 
using the MC simulation. 


The fraction of the photon momentum participating in the production of the two jets 
with highest is dehned 


as 


= 


1 


2pE, 


■{E 


,jetl -Jjje 


+ E 




( 1 ) 


where the variable y is given hy y = /s. The LO direct and resolved processes populate 

different regions in with the direct processes concentrated at high values of The 
variable was reconstructed via the above formula using the calorimetric-jet transverse 
energies and ycai = since many systematic uncertainties in the measurement of 

energy by the CAL cancel out event by event in the ratio of Eq. (|I]). The average resolution 
in was 0.05 for > 47 GeV and | cos 0*1 < 0.8. 


8 Acceptance corrections and systematic uncertain¬ 
ties 


The PYTHIA MC event samples of resolved and direct processes were used to compute 
the acceptance corrections to the dijet distributions. These correction factors also take 
into account the efficiency of the trigger, the selection criteria and the purity and efficiency 
of the jet reconstruction. The contributions from direct and resolved processes in the MC 
models were added according to a £t to the uncorrected x^^ distribution in the data. A 
good description of the Mil, | cos0*| and y data distributions was given both by PYTHIA 
and HERWIG. The differential dijet cross sections were obtained by applying bin-by-bin 
corrections to the measured distributions. The bin-by-bin correction factors differed from 
unity typically by less than 10 %. 

A detailed study of the sources of systematic uncertainty was performed. This study 
includes (a typical contribution for each item to the cross-section uncertainty is indicated 
in parentheses): 

• using the HERWIG generator to evaluate the correction factors for the observed dijet 
distributions (5%); 

• varying the cuts used to select the data while maintaining agreement between data 
and MG simulations ( 2 %); 
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• adding the contributions from direct and resolved processes according to the default 
cross sections as predicted by PYTHIA (2%); 

• using the LACl parameterisations of the photon PDFs for the PYTHIA MC samples 

( 2 %). 

The effects of uncertainties in the simulation of the trigger were negligible. All the above 
systematic uncertainties were added in quadrature to the statistical uncertainties. 

The absolute energy scale of the calorimetric jets in simulated events was varied by ±1%. 
The effect of this variation on the dijet cross sections was approximately =f 5%. This 
uncertainty is highly correlated between measurements in different bins. It is shown as 
a shaded band in Fig. [l|. In Figs. ^ to it has been added in quadrature to the other 
systematic uncertainties. An overall normalisation uncertainty of 1.6%, arising from the 
luminosity determination, is not included. 


9 High-mass dijet differential cross sections 

Using the selected data sample of dijet events, the differential dijet cross sections were 
measured in the kinematic region dehned by < 1 GeV^ and 134 < IF < 277 GeV. 
The dijet variables and cross sections were calculated using the two highest-!?;^* hadronic 
jets with > 14 GeV and —1 < < 2.5. For a given IF, events at high | cos0*| 

have smaller scattering angles and thus lower 77;^*. In order to study the |cos0*| and 
Ffjj distributions without bias from the 77;^* requirement, the cuts Mil > 47 GeV and 
I cos 0*1 < 0.8 were applied^. 

The cross-section da/dM^\ measured in the Mil range between 47 and 160 GeV and 
integrated over | cos0*| < 0.8, is presented in Fig. ^a). The data points are located at 
the mean of each bin. The measured da/dM^^ exhibits a steep fall-off over 3 orders 
of magnitude in the range considered. 

The cross-section da/dM^^ was also measured for the region > 0.75. The results 
are also presented in Fig. |I|(a) and exhibit a dependence on similar to that of the 
measurements integrated over the full range in The fraction of dijet events with 

^obs ^ increases with ilfjj from 57% at = 50 GeV to 98% at Mil = 139 GeV. 

To study the angular distribution as a function of the dijet mass, the normalised dijet 
cross section, \ja da/dx, was measured for 1 < x ^ 9 in four regions of M-ij; a is the 

^ The application of a lower cut on which is not close to the threshold avoids the infrared sensitivity 
[ pd| that otherwise would affect NLO QCD calculations for dijet production with symmetric cuts in 
the transverse energy of each of the two jets. 



dijet cross section integrated over the entire x range in each region. The normalised 
cross section has the advantage that the experimental and theoretical uncertainties are 
reduced, while allowing a precise test of the shape of the distribution in the calculations. 
The results integrated over the full range in are presented in Fig. and those in 
the region > 0.75 are shown in Fig. The measured l/cr da/dx decreases with 
increasing x in all the regions studied. 


9.1 Perturbative QCD predictions 


LO and NLO QCD calculations 0 are compared to the measurements in Figs. |l| to 
The calculations were performed using the GRV and CTEQ4 parameterisations 
of the photon and proton PDFs, respectively. The renormalisation and factorisation 
scales were chosen to be the highest as was calculated at the two-loop level using 
= 202 MeV, which corresponds to as(Mz) = 0.116. The calculations included only 
QCD hard-scattering processes and thus electroweak corrections such as qq —>■ qq via 
exchange or qq' —> qq' via exchange were not included. 


The uncertainty on the NLO calculations due to the absence of higher-order terms was 
estimated by varying the choice of the renormalisation and factorisation scales between 
/2 and 2E^^^ and amounts to less than 15% in the case of da/dM^K This uncertainty 
affects mainly the normalisation of the predictions and is therefore reduced for l/cx da/dx, 
for which it typically amounts to less than 5%. The uncertainty on the NLO calculations 
due to that on the gluon density of the proton was estimated by comparing with cal¬ 
culations based on the CTEQ4HJ parameterisations and amounts to less than 6% 
for da/dM^K The uncertainty due to the choice of photon PDFs was estimated by using 
the AFG parameterisations. The choice of parameterisation of the photon PDFs 
affects mainly the normalisation of the calculation. The normalisation for the case of 
AFG is smaller by approximately 10% than that based on GRV. All the above theoretical 
uncertainties were added in quadrature. 


NLO QGD calculations refer to jets of partons, whereas the measurements refer to jets of 
hadrons. An estimate of the effects of hadronisation was obtained by comparing the cross 
sections for jets of hadrons and jets of partons calculated with the PYTHIA program. 
The corrections for hadronisation effects were found to be within ±5% of unity both 
for da/dAfjj and 1/a da/dx', they therefore were neglected in the comparison of pQGD 
calculations to the measurements. 


9.2 Dijet mass distribution 

The LO and NLO QGD calculations describe the shape of the measured da/dM^^ well 
over the entire range of and see Fig. |^(a). The shape of this distribution in 
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the calculations is dictated by the x dependence of the photon and proton PDFs and 
by the dependence on the photon-parton or parton-parton centre-of-mass energy of the 
subprocess cross section. The LO QCD calculation of da/dM^^ is ~ 25% below the data. 
The inclusion of the NLO corrections signihcantly improves the description of the data. 
Although the data are still consistently above the calculations, the NLO QCD calculations 
are consistent with the data given the present theoretical (~ 15%) and experimental 
(~ 10%) uncertainties, as shown in Fig. 0(b). 

Since the region > 0.75 is dominated by direct processes, in which the photon behaves 
as a point-like particle, measurements restricted to that region in allow a test of the 
pQCD predictions with less influence from the photon PDFs. The NLO QCD predictions 
of dcr/dMjj have been compared to the data in the region > 0.75. In this restricted 
range, the NLO QCD calculations describe the data well both in shape and magnitude, 
as shown in Fig. |l|(c). 

Overall, no signihcant deviation between data and NLO QCD calculations is observed up 
to iVfjj = 139 GeV. 


9.3 Dijet angular distribution 


Figure shows that the shape of the NLO QCD calculation is in agreement with that of 
the measured l/cr da/dx in each region of dijet mass studied. An equally good description 
is obtained for the measurements restricted to x'^^ > 0.75, as can be seen in Fig. |^. The 
distribution in y reflects the spin of the exchanged particle: it is approximately uniform 
for two-body processes dominated by gluon exchange and proportional to 1/(1-l-y) for 
processes mediated by quark exchange. As an illustration, the LO QCD calculations for 
resolved and direct processes are shown in Fig. 0|: the x distribution for direct processes, 
which are mediated by quark exchange, is steeper than that of resolved processes, which 
are dominated by gluon exchange. The measured angular distribution is consistent with 
the pQCD description in terms of parton exchange up to the highest 


The ratio = ct(x^<x<9 ) useful to characterise, with a single number, the shape of the 
angular distribution in each range of dijet mass This ratio is shown in Fig. 0, with 
the value of xo chosen to be 4, so that is approximately 1. In this ratio, most of the 
systematic uncertainties cancel out, yielding a systematic uncertainty of less than 4% for 
all Mjj ranges considered. The uncertainties of the theoretical predictions for Ry. due to 
the choice of photon PDFs and from the choice of renormalisation and factorisation scales 
are also reduced. The NLO QCD calculations are consistently below, but compatible 
with, the data within the experimental (~ 8%) and theoretical (~ 10%) uncertainties. 
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9.4 Upper limits on Z®, and new heavy-resonance produc¬ 
tion 


The production of the electroweak gauge bosons and has been studied in the 
region 75 < < 100 GeV for the full ranged in The 1/a da/dx distribution 

is not consistent with a large contribution from and decays, as can be seen 
from Fig. |^, which shows that QCD processes are dominant. Simulations of both the 
QCD (non-resonant) dijet photoproduction background and the signal were used to 
hnd values that optimise the observation of a signal relative to the background for 
the upper cut on the measured value of | cos6**|, | cos6**|cut, and the dijet-mass window. 
The background expected from QCD processes was estimated by using the NLO QCD 
predictions with the normalisation obtained from a £t to the measured da/dM^^ in the 
region 47 < < 75 GeV. Since there is no evidence for a signal, an upper limit on the 

cross section for production has been derived. In the dijet-mass window of 91 ±9 GeV 
and for |cos0*| < 0.6, 230 events were observed while 223 events were expected from 
QCD processes. The acceptance for Z^ production after all the selection cuts was 14.6%. 
The resulting 95% C.L. upper limit in the kinematic region < 1 GeV^ is 


(^e+p^e+ZOX < 5.9 pb. 


In deriving the upper limit, the two different photon PDFs mentioned in Section 0 were 
used; the one that gave the most conservative limit was chosen. This upper limit on 
Z^ production is the first obtained at HERA. The Standard Model expectation for the 
process e~^p —> e'^Z^X integrated over all is 0.3 pb [^,Q. The same procedure was 
applied for production and an upper limit of ae+p^e+w^x < 7.4 pb at 95% C.L. 
was obtained. This limit is weaker than that already obtained by ZEUS using the 
leptonic decays The Standard Model expectation for the cross section of the process 
e+p —> e~^W^X is 0.95 pb • 


The increase of the yp centre-of-mass energy by an order of magnitude with respect to 
hxed-target experiments [Q allows a search for generic heavy resonances, denoted by V, 
with masses above 60 GeV in the reaction yp ^ P -|- X —jet -|- jet -|- X. While searches 
for such resonances formed by qq annihilation, such as Z', have been performed [Q with 
similar sensitivity in this mass range, there are models such as Technicolour which 
predict resonances^ produced preferentially by photon-gluon fusion. Such particles would 
not have been observed in earlier searches. 


A search for the photoproduction of a new heavy resonance decaying into two jets has 
been performed in the mass range 60 < M-p < 155 GeV. The method described above has 

Although the LO contribution to the production of and bosons is given by the resolved-photon 
processes (see Section ||), the resulting distribution peaks at ~ 0.75. 

® Some Technicolour models [Q predict the existence of light colour-octet, isospin-singlet pseudo- 
Goldstone bosons that would decay into two gluons. 
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also been used in this search to select the window for | cos0*| < 0.6. In this case, the 
background expected from QCD processes has been estimated by using the NLO QCD 
predictions with the normalisation obtained from a £t to the measured da/dM^^ in the 
region of values well below that of the signal. Since there is no evidence for a signal, an 
upper limit on the cross section times the branching ratio into two jets for V production, 
ae+p^e+px ■ Br('P —jet + jet), has been derived. The limit refers to the kinematic range 
defined by < i QeV^ 134 < IT < 277 GeV, > 14 GeV, -1 < V"* < 2.5 and 
Icos^*! < 0.6. The restriction to |cos0*| < 0.6 reduces the background from QGD- 
induced processes and the dependence on the spin and decay angular distribution of the 
heavy resonance (see Section 5). The upper limit on ae+p^e+px ■ By(V jet + jet) at 
95% G.L. is shown in Fig. ^ as a function of Mp. 


10 Summary 

Measurements of differential cross sections for dijet photoproduction have been made in 
e+p collisions at a centre-of-mass energy of 300 GeV using 42.7 pb“^ of data collected 
with the ZEUS detector at HERA. The dijet cross sections refer to jets identified with the 
kr cluster algorithm in the longitudinally invariant inclusive mode which were selected 
with > 14 GeV and —1 < < 2.5. The measurements were made in the kinematic 

region defined by < 1 GeV^ and 134 < IT < 277 GeV. The dijet differential cross 
section as a function of has been measured for > 47 GeV and | cos 0*1 < 0.8. 
Values of up to ~ 160 GeV are accessible with the present data. The NLO QGD 
calculations give a good description of the shape of the measured differential cross-section 
da/dM^k 

The dijet cross section as a function of y = (1 -l- | cos0*|)/(l — | cos0*|) has been measured 
in several Afjj ranges. The dependence of the distribution 1/a da/dx and of the ratio 
Rx = ff( 4 <x< 9 ) mass are found to be consistent with the predictions of NLO 

QGD in the range studied. The observed agreement in shape between data and the 
QGD calculations from Vfjj = 50 GeV up to iVfjj = 139 GeV conhrms the validity of 
the pQGD description of photon-parton interactions down to distances of approximately 
1.4 ■ 10-3 fm. 

The dijet events in the region 75 < < 100 GeV have been used to derive an upper 

limit on the cross section for production of ae+p^e+z^x < 5.9 pb at 95% G.L. for 
< 1 GeV^. Upper limits on the photoproduction of new heavy resonances decaying 
into two jets have been presented for 60 < Mp < 155 GeV. 
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Figure 1: (a) Dijet cross-section da/dM^^] the upper data points are for the 

full region in and the lower points are for the region > 0.75. The thick 
vertical bars represent the statistical uncertainties of the data, and the thin bars 
show the statistical and systematic uncertainties added in quadrature, except for 
that associated with the uncertainty in the absolute energy scale of the jets (shaded 
band). The LO (dashed line) and NLO (solid line) QCD parton-level calculations 
are shown. Both the measurements and calculations of da/dM^^, for > 0.75, 
have been multiplied by 0.1 for clarity of display; (b) the fractional difference 
between the measured da/dM^^ and the NLO QCD calculation integrated over the 
full region in (c) as (b), but for the region > 0.75. In (b) and (c), the 
hatched bands display the uncertainty of the NLO QCD calculations. 
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Figure 2: The normalised dijet cross sections, l/a da/dx, the full region of 
for four ranges in M^K The thick vertical bars represent the statistical uncertainties 
of the data while the thin bars show the statistical and systematic uncertainties 
added in quadrature. For comparison, NLO QCD parton-level calculations are 
shown as solid lines. The hatched bands display the uncertainty of the NLO QCD 
calculations. The LO QCD parton-level calculations for resolved (dashed lines) 
and direct (dotted lines) processes are also shown. The prediction of PYTHIA for 
the production of hF^ and bosons (dot-dashed line) is also shown in the region 
75 < M» < 100 GeV. 
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Figure 3: The normalised dijet cross sections l/cr da/dx, restricted to > 
0.75, for four ranges in Other details are as described in the caption to Fig. 
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Figure 4: The cross-section ratio, as a function of (a) inte¬ 

grated over the full region in (b) for the region > 0.75. The thick vertical 
bars represent the statistical uncertainties of the data while the thin bars show the 
statistical and systematic uncertainties added in quadrature. For comparison, the 
parton-level prediction of NLO QCD is shown as the solid line; the hatched band 
displays its uncertainty. 
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Figure 5: Upper limit at 95% C.L. on the cross section times branching ratio for 
the prodnction of new heavy resonances decaying into two jets, cxe+p^e+'px ■Br(P — 
jet + jet), as a fnnction of the resonance mass, M-p. 
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